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We review some known properties of the "emptiness formation probability" correlation which we have cal- 
culated numerically for spin-1/2 XX chains with constant (homogeneous) or alternating (dimerized) nearest- 
neighbor coupling and an external field (in z direction) for arbitrary temperature. The long-distance asymptotic 
behavior of this correlation is known to be Gaussian at zero temperature and exponential at finite temperature 
for the homogeneous chain. By simple analytical arguments the exponential behavior at finite temperature 
extends to the dimerized system. Numerical results for the dimerized chain confirm the exponential decay at 
finite temperature and show Gaussian decay at zero temperature. 
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1. Basic properties of the EFP correlation 

The "emptiness formation probability" (EFP) P(n) was first [ Q] defined as the probability 
of finding a "ferromagnetic island" of n sites in the ground state (or equilibrium state at finite 
temperature) of a spin-i antiferromagnetic chain. For obvious reasons it has also [[2] been called 
"ferromagnetic string formation probability" . From an experimental point of view P(n) is a com- 
plicated many-particle correlation. However, mathematically speaking, the EFP can be introduced 
as a natural elementary building block for constructing various correlation functions of integrable 
spin chains by using Bethe Ansatz techniques [ [3] , see also comments given in [0J. Depending on 
the situation at hand, a variety of analytical techniques was applied to derive more or less explicit 
results for the EFP or its long-distance asymptotics. The EFP can be represented as a multiple 
integral [OSj, as a Fredholm determinant [Q] or as a Toeplitz determinant [0[7]. Standard numer- 
ical techniques for one-dimensional quantum systems were used to obtain numerical data to which 
analytic or asymptotic results could be compared, for example the density matrix renormalization 
group (DMRG) at temperature T = [ [B] or quantum Monte Carlo simulations (QMC) at finite 
T[i. 

One model for which the EFP was discussed is the nearest-neighbor S = i XXZ antiferro- 
magnetic chain in an external field: 



H = jJ2 (s*S* +1 + SfSf +1 + A (sfSfa ~ D) -hj^Sf- (1) 



i=i 



For the special case A = (in which the model is known as the XX chain) a Gaussian decay of 
P(n) at long distances n was derived [[5] for T = 0, and it was conjectured that the behavior is 
Gaussian also for nonzero A. Subleading correction terms to the A = asymptotic behavior were 
found in [ [6] ; the present knowledge is summarized in the formula: 

n 2 / 1 \ 1 1 / 1 \ 1 

lnP(n) = T lnf— (J + h)\ -- Inn --In I— (J- h)j +^2 + 3^-1) + o(l). (2) 
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(Here £(x) denotes the Riemann zeta function.) Other asymptotic formulae were derived [[2] for 
h values close to the saturation fields h — ± J, where the ground state of the XX chain becomes 
a completely polarized ferromagnet. A related T = correlation function involving a weakly 
ferromagnetic string of spins was calculated [ [9] using bosonization methods, and the value of 
P(6) was calculated [\W\ for the XXX model (A = 1) at T = h = 0. For arbitrary |A| < 1, an 
asymptotic formula (for T — h — 0) of the form 

P(n) Pa An-^G^ (3) 

with explicit expressions for C and 7 (as functions of A) was suggested [111]. This formula repro- 
duces the known exact results for A = 0, | and compares favorably to numerical data for other A 
values. 

For nonzero temperature the EFP correlation of the XXZ chain is expected to show exponential 
decay for sufficiently large n. The following simple argument for this asymptotic form was given 
by Boos and Korepin [[12]. The EFP can be defined as the expectation value of the projection 
operator onto a state with n adjacent up spins: 

/" \ Tre~& ft" Pj 



with 



Pi = S' j+ \. (5) 



The partition function for an iV-site chain is 



Z N = Tre~—=e~— (6) 

in the limit of large N, where / is the free energy per site. Note that the saturated ferromagnetic 
state with all spins up is a zero-energy eigenstate of the Hamiltonian ([I]). (Let us temporarily 
consider the case h = for simplicity.) The product of projection operators in the numerator of 
P(n) (TJ| thus projects onto states with n up spins, for which the corresponding n-site "partial 
Hamiltonian" can be replaced by zero. In other words, neglecting boundary effects, the numerator 
of P(n) is nothing but the partition function of an (N — n)-site chain: 



The neglect of boundary effects in the above "derivation" becomes increasingly problematic at low 
T, because the zero-temperature correlations in the model ([T]) have long-range power law decays. 

A finite magnetic field h does not change the above argument essentially; of course the free 
energy per particle depends on h. The ferromagnetic state still is an eigenstate of the Hamiltonian, 
but its energy eigenvalue changes from zero to —h/2 per site; accordingly the asymptotic form of 
P(n) is 

/ n ( „,„ , „ h 



P(n)~exp\ — \f(T,h) + ~\). (8) 

It should be noted that nowhere in the above line of argument we exploited the fact that the 
nearest-neighbor exchange interactions J in the Hamiltonian ([T]) are all equal; the asymptotic 
formula |(5J| thus should hold also for modulated (for example alternating) exchange interactions. 
However, for modulations with very long wavelength or for random couplings the concept of a 
constant free energy per site loses its meaning. 

In the XXZ chain ^ the total spin z component S* ot = Yli=i &t i s a conserved quantity. That 
is different for the (anisotropic) XY chain [ [T31 [Tl] 

N , v N 

H = J E ( (! + 7)S?5^i + (1 ~ 7)^5?+! ) - h £ SI (9) 
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with anisotropy parameter 7. The (7, /i)-plane shows several ground-state phases. The XX model 
discussed above corresponds to the line 7 = 0, where the ground-state EFP is asymptotically 
Gaussian [ 6 . Away from that line, however, the EFP is asymptotically exponential [ [T5l [TBI [7], 
For nonzero temperature the EFP is asymptotically exponential for arbitrary 7 and h [ 16J. 

A model which shows surprising similarities, but also important differences to the anisotropic 
XY chain © is the dimerized XX chain 

N , v N 

H = Jj2(l-(-lYt) lS?S* +1 +SVS? +1 )-hY,S?. (10) 

i=l ^ ' z=l 

The parameter S quantifies the degree of dimerization; for S = 1 the system decomposes into 
decoupled spin dimers. Ground-state properties (in particular, dynamic correlations) were studied 
in parallel for both models, © and (fT0|) . by Taylor and Miiller [[17]. Dynamic properties at finite 
temperatures were studied in [ I18j for the dimerized model and in [ I19j for the anisotropic model 
(including an additional Dzyaloshinskii-Moriya interaction). Both models map to noninteracting 
Fermi quasiparticles under the Jordan- Wigner transformation; see next section for details. For 
7 = in ([9]) and S = in (jXOj) the two models are identical, with a gapless single-particle energy 
spectrum for the Jordan- Wigner fermions. For small nonzero 7 or S and moderate magnetic field h a 
gap appears in the single-particle spectra of both models. However, the parameter dependence of the 
energy spectra and the physical meaning of the quasiparticles are fundamentally different between 
the two models. In the dimerized model, the number of quasiparticles is a conserved quantity strictly 
related to the total spin z component Stot- In the anisotropic model, in contrast, Stot is no longer 
conserved and the Bogoliubov rotation involved in the diagonalization of the model mixes creation 
and annihilation operators (corresponding to spin raising and lowering operators). Consequently 
the ground state of this model is very simple in terms of the Jordan- Wigner quasiparticles, while it 
is a complicated superposition involving many Stot eigenstates in terms of the original spins. The 
shape of the single-particle spectrum in the dimerized model depends only on S (J fixing the overall 
energy scale). The spectral gap separates two bands containing equal numbers of single-particle 
states. The magnetic field h (equivalent to a chemical potential) determines the occupations of the 
bands, or, in spin language, the total magnetization Stot- In the anisotropic model, in contrast, 
the shape of the single-particle spectrum depends on both h and 7, the two energy bands contain 
different numbers of states, the lower quasiparticle band is always filled, and the upper band is 
always empty, see [ 117] for more details about the ground states of the two models. 

It is thus of interest to study the EFP correlation for the dimerized model and to compare the 
results to those obtained [Q~51[TH1[7] for the anisotropic model. Our results show that the EFP of the 
dimerized system is asymptotically Gaussian at T = and exponential at T > 0. The exponential 
behavior is described by the simple free energy argument leading to (|SJ). 

The remainder of the paper is organized as follows. In section [2] we explain the numerical pro- 
cedure which can be used to calculate P(n) for arbitrary inhomogeneous XX chains. Section [3] 
contains numerical results for the homogeneous chain which are compared to the expected asymp- 
totic behavior. Finite-size and boundary effects are also discussed. Section @] presents results for 
the dimerized chain, with two different values for the dimerization parameter. 



2. The numerical procedure 

We consider the general S — 1/2 XX chain 

jV-l N 

h=J2 Ji(SfSf +1 + SfSf +1 ) -hJ2s? (11) 

i=l i=l 

which is one of the simplest quantum many-body systems conceivable, because many of its prop- 
erties can be derived from those of noninteracting lattice fermions. 
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The Hamiltonian (fTT|) describing an open-ended TV-site spin- 1/2 XX chain can be mapped to 
a Hamiltonian of noninteracting fermions, 

N-l N 
H f = 5 E J ^+l + C Ul*) ^(cU - -) (12) 

i=l i=l 

by means of the Jordan- Wigner transformation [I13| 120] between spin and Fermi operators: 

Pi = St + \ = 4a, (13) 

St = C K] = H(l - 2 c t Ci )4. (14) 

(=1 

In the homogeneous case Ji = J, the fermion Hamiltonian is 

^ = ^^4^ + — (is) 

fc 

where the operators c^, and Ck create and destroy a fermion in a one-particle eigenstate, respectively. 
The one-particle energy eigenvalues are 

VK 

£h = J cos k — h, fc = ,v = l,---,N (16) 

N + 1 

and the eigenvectors are sinusoidal functions of the site index i. For a dimerized chain, where 
Ji alternates between even and odd i (compare (I10|) ). the eigenvalues £/- and the corresponding 
eigenvectors are known analytically [ 17J. In the notation of the Hamiltonian (flQ|) the one-particle 
energies are given by 

Ek = —h + Jsign(cos k) V cos 2 fc + S 2 sin 2 fc. (17) 

It should be noted that (fl~7|) was derived in the thermodynamic limit N — > oo. For finite it 
is important to distinguish between even and odd chain lengths in the dimerized case. Detailed 
discussions of both cases can be found in [ [211 E2]- The data presented in the present paper 
were obtained by numerical diagonalization without making explicit use of the dispersion relations 
(fT6TT7l) . 

The ground state has all single-particle states with negative energies occupied by Jordan- Wigner 
fermions while all other states are empty. For \h\ < J 5 the zero-energy level lies within the spectral 
gap of (|17| and hence the ground state does not depend on h in this range. The same is true 
for \h\ > J where the ground state is either completely occupied by Jordan- Wigner fermions or 
completely empty In the intermediate field range, J 5 < \h\ < J, the groound state contains a 
partially filled band of Jordan- Wigner fermions. Consequently we may expect behavior similar to 
that of the homogeneous chain for \h\ < J. 

For a trimerized system, that is, Ji+3 = J,: in (llip the eigenvalues and eigenvectors still are 
known analytically [[53], whereas for periodically varying couplings with a larger period p they can 
in general only be obtained by numerically solving apxp eigenvalue problem; see, for example 
[ 124] . For general Ji neither eigenvalues nor eigenvectors are available analytically, however, both 
are easily obtained from the solution of a tridiagonal eigenvalue problem with standard numerical 
procedures [125]. 

The asymptotic formula Q for the EFP correlation may be evaluated more explicitly. Taking 
into account (|15[) and the standard formula for the free energy of free fermions we obtain 



P(n) = c(T, h) exp -— — V ln(l + e"^) - - . (18) 




This approximate asymptotic formula for the EFP correlation should hold as long as the concept 
of a spatially (roughly) constant free energy per site makes sense, as discussed above. We have 
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introduced a temperature- and field-dependent prefactor c(T, h). For the homogeneous XX chain, 
Shiroishi et al. [ [3] have derived an explicit expression for c(T, h) . 

The exact numerical evaluation of the EFP correlation proceeds as follows. Due to (fT3|) 

P(n) = (cjci4c2 • • • C^Cn). (19) 

This many-fermion expectation value can be evaluated using the Wick-Bloch-De Dominicis Theo- 
rem [ 26 a . Mathematically this means that the expectation value (|T9|) can be expressed as a Pfaffian 
involving only elementary expectation values. Pfafhans are close relatives of determinants and we 
refer the reader to the literature for their properties [ [STJ . The numerical evaluation of Pfafhans 
proceeds along similar lines as that of determinants. Elements can be reduced to zero by operations 
which are known to leave the value of the Pfaffian invariant. After production of sufficiently many 
zero elements the evaluation of the Pfaffian becomes trivial due to an expansion theorem. An im- 
plementation along these lines was described by Derzhko and Krokhmalskii [[28]. We use a similar 
algorithm here. A recursive scheme for evaluating Pfafhans was used by Jia and Chakravarty [ 29J. 



3. Results for the homogeneous case 
3.1. T — 




We have tested our numerical evaluation for chains with N = 256. Figure [T] shows the EFP 
correlation P(n) for a string of spins starting at site i = 110, for several values of the magnetic 
field h. 

In order to assess the influence of finite-size and boundary effects we have also performed 
computations for i — 111, and also for N = 512 and i = 221. We found that at T = there are still 
finite-size and even/odd effects of the order of up to one percent. Given the fact that ground-state 
correlations in quantum spin chains tend to be long-ranged, this does not come as a surprise. 
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On the scale of the figure there is near perfect agreement between numerical evaluation of 
the correlation and the asymptotic formula as already noted by Shiroishi ct al. [ 6 . Only 
at extremely small values of the correlation function (outside the range of the figure) there are 
significant deviations from the asymptotic behavior. These are, however due to precision problems 
in the numerical evaluation of P(n). 

3.2. T > 




Figure 2. The EFP correlation of the homogeneous chain at T = 0.1, for magnetic fields h — 0.99, 
0.5, 0.2, 0.1, 0, -0.1, -0.2, -0.5, -0.99 (top to bottom). Crosses represent the numerical evaluation, 
solid lines show the asymptotic behavior ((18}, where we have put c(T, h) = 1. 

We have numerically evaluated the EFP correlation P(n) for N = 256 chains for a string of 
spins starting at site i — 90 for several values of T and the magnetic field h. Figures 1213141 and [5] 
show results for T ~ 0.1, 0.5, 1, and 10, respectively. 

In all cases it was found that the asymptotic formula (fTS)) describes the numerical data very 
well. The prefactor c(T, h)(> 1) is large at low T and negative values of h, as already observed by 
Shiroishi et al. [[6] (compare their Figure 6). For intermediate and high temperatures we observed 
a very small but seemingly systematic slope when we plotted the difference between the logarithm 
of the numerical value of P(n) and the expression in the exponent of (fT8|) versus n. Due to the 
smallness of the effect we did not further pursue the origin of this difference, however. 

We have also assessed the influence of finite-size and boundary effects on P(n) for T — 0.1 
and h = ±0.99, ±0.1, and 0. For these parameter values we compared the numerical results for 
N = 256, starting sites z = 90 and 91, and N = 512, i — 240. No differences beyond the 10 -7 level 
were found. 

In the limit T — ► oo P(n) = (i) for arbitrary h. That is precisely what the h = results 
already show at T = 10 (see Fig. [5]). For some selected h and T values we have also compared our 
values for the decay length of P(n) to those given in Figure 5 of Shiroishi et al. [[B]. The results 
agree. 

At low temperature a crossover from the Gaussian T = behavior ([2|) at small n to the 
exponential asymptotics (I18|) at larger n can be observed, as shown in Fig. [6] 
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Figure 5. Same as Figure [2] for T = 10. 




n 



Figure 6. The EFP correlation of the homogeneous chain at T — 0.05, for magnetic fields 
h — 0.99, 0.5, and 0.2 (top to bottom). Crosses represent the numerical evaluation, solid (red) 
lines show the exponential asymptotic behavior (|18[) , with adjusted value for c(T,h), dashed 
(blue) lines show the Gaussian asymptotic behavior ((2} for h — 0.5 and 0.2. 
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n 

Figure 7. The EFP correlation P(n) of the dimerized chain as a function of n 2 for T — 0, S = 0.2 
and h = -0.75, -0.5, -0.25, -0.24, -0.23, -0.22, -0.21, -0.20, +0.201, 0.202, 0.203, 0.204, 0.205, 0.21, 
0.25, 0.5, and 0.75 (bottom to top). Heavy lines denote the data for h = -0.75, -0.5, -0.25, -0.20 
(dashed), +0.205, 0.25, 0.5, and 0.75. On the scale of this figure the data for h = —0.20,0, and 
+0.20 coincide. Note the rapid variation for h slightly larger than 0.20. (The data for h = 0.202 
and 0.203 coincide on this scale, as do the data for h — 0.204 and 0.205.) 



4. Results for the dimerized case 

We have calculated P(n) for weak and strong dimerization, 5 = 0.2 and S = 0.8. 
4.1. S - 0.2, T — 

For given dimerization r5 > there are two critical values for the magnetic field, h = ±<5. For 
\h\ < 5 the ground state does not depend on h. In the Fermion picture it consists of a completely 
filled lower band and a completely empty upper band. It may thus be expected that correlation 
functions do not change with h for h within this range. 

Figure [7] shows P(n) for a large number of h values, with special attention to the critical 
regions near h — ±5. The plot of lnP(n) vs n 2 convincingly demonstrates an overall Gaussian 
behavior, but a closer look at the region n < 10, for example, reveals some oscillations due to the 
dimerization. As expected, the EFP correlation does not change as h varies between —S and S. 

Figure [8] shows the values of the Gaussian decay parameter a 2 obtained from the data of Figure 
[7]by a least-squares fit of the form In P(n) = ao + a\n — a^r? . The lines between the data points are 
guides to the eye only. The band-edge singlarity at h — S is very prominently visible. In contrast, 
the behavior near h = —5, where the lower quasiparticle band approaches complete filling, is less 
singular. 

The changes caused by varying the initial site i are shown in Figure [9] for a sub-critical field 
h = 0.2 and a super-critical field h = 0.5. The solid and dashed lines are for i = 110 and 111, 
respectively. The P(n) values for even and odd i seem to coincide for odd string lengths n. A closer 
look at the data shows that at h = 0.2 this coincidence of data for odd n is indeed perfect. For 
h = 0.5 differences at the 1% level show up for all odd n. This behavior is not surprising since for 
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Figure 8. Gaussian decay parameter 0,2 obtained from the numerical P(n) data shown in Fig. [7] 
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Figure 9. The EFP correlation P(n) of the dimerized chain for T = 0, S = 0.2, h — 0.5 (upper 
set of curves) and h = 0.2 (lower set of curves). The strings of sites considered start at site 
i — 110 (solid lines) or i = 111 (dashed lines) of a N = 256 chain. 
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Figure 10. The EFP correlation P(n) of the dimerized chain for T = 0.01, 8 = 0.2, and h = 1.2, 
1, 0.5, 0.2, 0.1, 0, -0.1, -0.2, -0.5, 1-, -1.2 (top to bottom). Symbols correspond to numerically 
calculated values, solid lines to the asymptotic formula 1181 On the scale of the figure the data 
for h = 1 and 1.2 cannot be distinguished. 



odd n the string of n sites contains equal numbers of "weak" and "strong" bonds for both even 
and odd i, hence the EFP as a global quantity of the whole string should be the same, at least in 
the limit of infinite system size. For even n, however, differences arc to be expected, and arc indeed 
present in the numerical data. As an extreme example, consider n = 2: P(2) refers to a single bond 
and will significantly depend on whether that bond is strong or weak. 

In order to assess the effects of system size we have also compared the N = 256 data to N = 512 
data (for i = 220 and 221). The data for h = 0.2 do not depend on N, whereas the h = 0.5 data 
vary at the percent level between the two N values. 

As we are dealing with finite open systems, boundary effects are an issue to be considered. 
For a given absolute value of the dimerization parameter \S\ the chain may either start with a 
"strong" bond or with a "weak" bond. It turns out that switching back and forth between these 
two possibilities (at N — 256) is equivalent to switching between even and odd i- The equivalence 
is in fact perfect for subcritical values of \h\, for example h — and 0.2. At h — ±0.5 we observe 
, differences at the percent level, the size of which does not change significantly between N = 256 
and N = 512. 

It should be an interesting task to generalize the asymptotic formula J5} to the dimerized case 
along the lines of [ [5] in order to better understand the behavior observed in Fig. [5] 

4.2. S = 0.2, T > 

Figure [10] shows data for the low temperature T = 0.01 and a wide range of h values. The 
strings of sites considered again start at site % = 110 of a N = 256 chain. Comparing the numerically 
calculated values of P(n) to the asymptotic formula (Tig]) we see clearly that at this low temperature 
the EFP correlation for n < 20 still displays the typical Gaussian zero-temperature behavior. The 
correlations for h = —0.2, -0.1, and stay together very closely for n < 20, whereas those for 
h = 0.1 start to deviate for n > 10 and those for h = 0.2 even earlier. This is in contrast to T = 
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n 



Figure 13. Similar to Fig. [TT] for T = 0.3 

(see Fig. [7]) where the correlations for all \h\ < 0.2 coincide because the ground state does not 
change in this h range. 

This behavior changes as T increases, as shown in Figures [TT] , [T2l and [I3j for T — 0.03, 
0.1, and 0.3, respectively. Figure [TT] shows a crossover of P(n) from Gaussian behavior at small 
n to exponential behavior at larger n. However, at this rather low T value the prefactor c(T, h) 
in the asymptotic formula (|18|) is still quite large. Superimposed on the asymptotic behavior the 
numerical data at T = 0.03 and h = 0.5 show the modulation expected due to the alternating 
exchange couplings. That modulation becomes more visible and the prefactor c(T, h) approaches 
unity as T grows, see Figs. Upland [TBI 



4.3. S - 0.8 

For strong dimerization we present data for h = ±0.99, ±0.9, ±0.8, and 0. Figure [Ml shows the 
results for T = 0. The plot of In P(n) as a function of n 2 shows that the decay is still Gaussian apart 
from the expected modulation due to the dimerization, which for <5 = 0.8 is of course expected to 
be more visible than for <5 = 0.2. Note that the data points for h = 0, ±0.8 coincide as the ground 
state stays the same for thes h values. 

At T = 0.03 (Figure [T5]) the data seem to fall naturally into three groups, namely h < 0, h = 0, 
and h > 0. The h < data do not seem to reach the proposed asymptotic behavior ([T8| within 
the range n < 20, whereas the other data do so quite nicely. This is not true, however, as closer 
inspection of the numbers reveals. The data for h < reach the asymptotics later and with a much 
larger prefactor c(T, h) in (fTS]) . 

At T = 0.1 (Figure [TB]) all n > 5 data fit to the expected exponential behavior apart from the 
expected modulation due to dimerization. Note that at this 'high' temperature the h — data 
stand out alone and the data for the critical field values h — ±0.8 are much closer to the data for 
h = ±0.9 and ±0.99, in contrast to the situation at T = (Fig. [HI). 
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Figure 14. The EFP correlation P(n) of the dimerized chain for T = 0, S = 0.8, and h = 0.99, 
0.9, 0.8, 0, -0.8, -0.9 and -0.99 (top to bottom). The data for h = 0, ±0.8 coincide. 
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Figure 15. The EFP correlation P(ri) of the dimerized chain for T = 0.03, 8 = 0.8, and h = 0.99, 
0.9, 0.8, 0, -0.8, -0.9 and -0.99 (top to bottom). Symbols denote numerically calculated values, 
solid lines show the asymptotic formula (JTSJ, with c(T, K) = 1. 
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5. Concluding Remarks 

We have examined the emptiness formation probability (EFP) correlation for homogeneous 
and dimerized spin-1/2 XX chains numerically by evaluating Pfaffian forms. For the homogeneous 
XX system the EFP is known to show asymptotically Gaussian behavior in the ground state 
and exponential behavior at finite temperature. For a homogeneous (anisotropic) XY chain the 
EFP is known to be exponential at both zero and finite temperatures. The anisotropic XY chain 
has a gapped single-particle energy spectrum, as does the dimerized XX chain. However, the 
EFP correlations of the two chains behave differently. Our calculations show that the EFP of the 
dimerized chain is Gaussian at at zero temperature and exponential at finite temperature. The 
exponential behavior at finite temperature can be interpreted within a simple free energy picture 
developed for the homogeneous chain. The generalization of the T = asymptotic expansion from 
the homogeneous to the dimerized system is left as a task for future work. 
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